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1. Introduction 



Neutrinoless double beta decay {0^(3(3) is a lepton number violating process, which, if observed, 
would prove that neutrinos are Majorana particles [1]. New physics beyond the standard 
model is required to make the process observable, and there are several different theoretical 
frameworks that could provide the necessary operators (see the review in Ref. ['^j). One of 
those theories is the left-right symmetric model (LRSM) [3-~], in which parity is restored at 
high energies and right-handed neutrinos are naturally included as part of an SU{2) doublet 
of the extended gauge symmetry. In that case there are a number of new physics contributions 
to Oz//3/3, either from right-handed neutrinos or Higgs triplets, with the rate for double beta 
decay linked to neutrino mass. This connection can be both indirect, through the couplings to 
and/or mixing with right-handed neutrinos, as well as direct, via the standard light neutrino 
contribution (see Ref. [8] for one of the first discussions of Oz//3/3 in the LRSM). 

In the simplest version of the LRSM one expects the scale of parity restoration to be rather 
high, i.e., around the GUT scale of 10^^ GeV. Indeed, if all couplings in the scalar potential of 
the theory are of order one then this conclusion follows naturally [7]. Nevertheless, there is still 
enough freedom in parameter space to allow one to consider TeV-scale left-right symmetry, 
which leads to several distinct and observable signatures in present-day experiments probing 
leptonic processes. On the other hand, the quark sector of the TeV-scale model is severely 
constrained, due to the presence of flavour changing neutral currents (FCNCs) induced by the 
neutral components of Higgs bi-doublets that are introduced to break electroweak symmetry. 
These affect meson mixing, CP violation in meson decay and the neutron electric dipole 
moment, and one needs the neutral component of the Higgs bi-doublet to be heavier than 
about 25 TeV to avoid conflict with experiment. The mass of the right-handed VT-boson 
{Wr) can however still be around 3 TeV, and current LHC data is already beginning to probe 
Wpi masses of this order [9, 10]. Indeed, the latest limits from the CMS experiment are roughly 
~ 2.5 TeV (see Fig. 6). With right-handed neutrinos of similar mass or lighter there 
are observable effects in OufSP and lepton flavour violation (LEV). The connection between 
double beta decay, LHC and lepton flavour violation has recently been studied by several 
authors [11-17]. 

From the theoretical point of view, the LRSM provides a natural framework for both the 
type I [18-22] and type II [6, 23-27] seesaw mechanisms, mediated by right-handed neutrinos 
and Higgs triplets, respectively. In this way the smallness of neutrino mass is connected to the 
restoration of parity at high energies, and the OuPP process can proceed via the same mediators 
that lead to neutrino mass. It is however rather difficult to pin down the mechanism by which 
the process occurs. A simplified case that has already been studied in the literature is that 
of type II seesaw dominance for rrii, [11], which restricts the number of parameters by making 
the right- and left-handed Majorana mass terms proportional to each other. We perform a 
detailed investigation of this case and include LEV constraints explicitly in the calculation of 
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the 0u/3(3 half-life. 

The case of type I seesaw dominance is more complicated: there are some contributions to 
Oz//3/3 that involve the left- and right-handed sectors individually as well as others that involve 
both sectors, through "left-right mixing". A simplified version was studied in Ref. [1 J], and a 
useful formula relating the various mass matrices of the theory was presented in Ref. [1 ')], for 
the case of symmetric Dirac coupling. Since the left-right mixing is always a ratio of the Dirac 
and Majorana mass scales, Oz//3/3 processes involving left-right mixing can be enhanced for 
specific Dirac mass matrices. This enhancement [28, 29] is also required for collider signatures 
of the TeV-scale seesaw mechanism (see the review in Ref. [-jU]), and there have been several 
studies of related phenomenology [31-34]. In the LRSM case both the so-called A- and i]- 
diagrams could give large contributions, although the latter is further suppressed by the mixing 
between left- and right-handed gauge bosons. This idea has also been discussed in the context 
of the inverse process e~e~ — )■ Wj^W^ [35] and was further emphasized in extended seesaw 
versions of the LRSM [16, 17]. We perform a thorough analysis of the type I seesaw scenario, 
paying attention to the correct nuclear matrix elements for the different diagrams as well as 
the often severe constraints from lepton flavour violating phenomena. 

The paper is outlined as follows: in Section 2 we briefly summarize the theoretical details of 
the left-right symmetric model (the reader familiar with the LR model may skip this section), 
and in Section 3 we provide a detailed discussion of the Oi//3/3 and LFV processes in the model. 
Section 4 is a quantitative analysis of the various Oz//3/3 amplitudes in the limit of type I or 
type II seesaw dominance; we summarize and conclude in Section 5. Details of decay widths 
and loop functions for LFV processes can be found in the appendix, which the reader may 
skip as well. 



In the left-right symmetric model, the Standard Model is extended to include the gauge 
group SU{2)r (with gauge coupling qr ^ ql), and right-handed fermions are grouped into 
doublets under this group. Thus we have the following fermion particle content under SU (2)^ x 



2. The left-right symmetric model 



SU{2)r X U{1) B-L- 




(1) 





with the electric charge given by Q = T| + T| + and i = 1,2, 3. The subscripts L and R 
are associated with the projection Pl,r = |(1 =F 75)- In order to break the gauge symmetry 



3 



and allow Majorana mass terms for neutrinos one introduces the Higgs triplets 

with Ai ~ (3, 1, 2) and ~ (1, 3, 2); the electroweak symmetry is broken by the bi-doublet 
scalar 

(2,2,0). (4) 




The relevant Lagrangian in the lepton sector is 

C^Y^- L^L{f(t> + f^)L'R - tlia^/^LhLL'L - L%ia2ARhRL'ji + h.c, (5) 

where = (T2(j)*(T2'i f,g and Hl^r are matrices of Yukawa couplings and charge conjugation is 
defined as 

= C^^Ir = {r)R,L , C = Z7270 ■ (6) 

If one assumes a discrete LR symmetry in addition to the additional gauge symmetry, the gauge 
couplings become equal {gi = gR = g) and one obtains relations between the Yukawa coupling 
matrices in the model. With a discrete parity symmetry (L^ -H- Lr, (f) (p^ -H- A^) it 
follows that = h*j^, f = P , f = P\ with a charge conjugation symmetry (L^ -H- {Lr)'^, 
(f) -H- 0^, Al <h- Ar) we have h = hi = hR, f = f^, f = f^. Applying these symmetries 
simplifies various expressions in the model, as will be discussed later. 

Making use of the gauge symmetry to eliminate complex phases, the most general vacuum 

is 

After spontaneous symmetry breaking, the mass term for the charged leptons is 

>Cl.a.s = -4W^ + h.C., (8) 

where the mass matrix 

Me^^(K2e'^f + Kj) (9) 

can be diagonalized by the bi-unitary transformation 

4,« = VlReL,R , vPMeVl = diag(me, m^, m,) . (10) 

With a discrete parity (charge conjugation) symmetry, becomes hcrmitian (symmetric), 
so that the condition ¥[ = (Vl = V^*) holds. In the neutrino sector we have a type I + II 
seesaw scenario. 



Cass = -XMX + h.c. = -i (u', z.-) I I ^ I +h.c., (11) 
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with 

= i=(/€i/ + «:2e-*"/) , ML = V2vLe''^hL, MR = V2vRhn. (12) 

Again, with a parity (charge conjugation) symmetry we have = Mj^ (M^ — ^d)- ^^le 
most general case the phase 6l cannot be set to zero, but in the type II dominance case we will 
study it is simply an overall phase and has no effect on the resulting neutrino mass matrix (in 
the type I dominance case it plays no role since vl = 0). Due to the presence of the so-called 
"VEV seesaw" relation relating the various VEVs, one expects x = vlVr/k\ = 0{1), since 
a; is a function of (order one) couplings in the scalar potential [7]. However, from a purely 
phenomenological point of view, x can take any value between and 10^^ [ ]. Assuming 
that Ml <C <^ Mr, the light neutrino mass matrix can be written in terms of the model 
parameters as 

m, = ML- MdM^^MI = V2vLe''^hL - ^^hohr^^hl , (13) 

where 

hD = —^ , = Kl + K2 . (14) 

The symmetric 6x6 neutrino mass matrix M^ in Eq. (11) is diagonalized by the unitary 6x6 
matrix [37-39] 

W) \T V) \ -fit l-\B}R)\ti Vn) 

to W^MyW* = diag(mi, m2, m3. Ml, M2, M3), where the unitary matrices and Vr, are 
defined by 

Ml - MdM^^MI = V, diag(mi, m^, mg) Vj, 

1 16) 

Mr = Vr diag(Mi, M2, M3)KJ , 

and the matrix R = Mj^M^^ + 0{M%{M^^f) describes the left-right mixing. The neutrino 
mass eigenstates n = til + n\ = rf are defined by 




n'r = ^. = Wn, 



Note that the unitarity of W leads to the useful relations 

YlYl^ = UU^ + SS^ = 1 = V^V^^ = TT^ + VV^ and V}^V^^ = UT^ + S0 = 
with the unitary 3x6 matrices = {U S) and = (T V) defined in Eq. (15). 



;i7) 
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The leptonic charged current interaction in the flavour basis is 



where 



^% = ^ W^'Pli^'^Z, + ^'i^Pri^'W^,] + h.c, (19) 





(20) 



cos ^ sin ^ e 
sin ^ e~*" cos ^ 

characterizes the mixing between left- and right-handed gauge bosons, with tan 2^ = — jTi^l • 
With negligible mixing the gauge boson masses become 

- ^ , and ^ ^ ' (^1) 

and assuming that^ ^2 < ^i, it follows that 

{==-W4--2^(^)'. (22) 

so that the mixing angle ^ is at most'^ the square of the ratio of left and right scales {L/Ry. 
Here we assume L ~ 10^ GeV corresponds to the electroweak scale and R ~ TeV to the scale 
of parity restoration, vr- For small ^ the charged current in the mass basis becomes 

C'^l, = -j= [hrKLULiWr, + ^e^^W,-^) + T^YKRuli-^e-^^Wr, + W,-^)] + h.c. (23) 

Here and Kji are 3x6 mixing matrices 

Kl = vIWZ, and KR = V^Wli*, (24) 

connecting the three charged lepton mass eigenstates ii to the six neutrino mass eigenstates 
(ui^Nif, {i = 1,2,3), with [using Eq. (18)] KlKI = KrK^^ = 1 and KlK^ = 0. The 
standard neutrino mixing matrix is just the left half of Kl, i.e., f/pMNS = ^lU. 

In this model one also expects a new neutral gauge boson, Z', which mixes with the standard 
model Z boson. The mass eigenstates Zi^2 have the masses 



9 , gcosOw /2cos2 6'vy 

^z, ^ 7^ ~ — , and mz^ ~ ^==vr ~ W — — mw^ , 25 

2cosfc'vK cosdw vcos26'h^ V cos26'vi/ 



where 9 = e/ sinOw and the U{1) coupling constant is g' = e/\/cos29w Again one expects 
the mixing to be of order (L/R)'^, i.e.. 



1 . _i g'^f^lVcos 29w Vcos 29w / T^f^z^ 



2 



2 2c^(m|^-m|J - "^Zi V^^Za/ 

Eqs. (21) and (25) imply that — l.Tmiy^. The current limits [41, 44] on the neutral gauge 
boson parameters are niz' > 1.162 TeV and |0| < 1.2 x 10~^. In addition, the current limits 
on the doubly charged triplet masses are [15] 'mg±± > 409 GeV and mg±± > 322 GeV. The 
theory predicts m^±± ~ vr, assuming order one couphng constants in the scalar potential. 



^This is justified if one assumes no cancellations in generating quark masses [40]. 

^Although the experimental limit is ^ < 10^^ [ ], for tovKr — O(TeV) one has ^ ^ 10""^ [42]; supernova 
bounds for right-handed neutrinos lighter than 1 MeV are even more stringent < 3 x 10"^) [42, 43]. 
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3. Oz//3/3, lepton flavor violation and collider physics 



3.1. Neutrinoless double beta decay 
3.1.1. Particle physics amplitudes 

Here we summarize the various possible diagrams for Oz//3/3 in left-right symmetric models (for 
one of the first analyses on this topic, see Ref. [S]). The Lagrangian in Eq. (23) can be written 
as 



^ i=l 



^ i=l 

+ e7?7'^(T,>£, + V:,Nn^{-ie-''^W{^ + W^^\ + h.c. 



(27) 



where in the second line we have assumed a basis where the charged leptons are diagonal (we 
will use this basis from now on, thus expressing all processes in terms of the matrices U , S, 
T and V). Oz//3/3 amplitudes arise from second order terms in perturbation theory: it is clear 
that one can combine either two left-handed currents, two right-handed currents or one left- 
and one right-handed current. The relevant mixing matrix element also depends on whether 
light or heavy neutrinos are exchanged in the process; the matrices U, V, S and T are (to 
second order in Mij/Mji) 



U 



l--M^M^\M^Mj 



K, V 



Vr, 



{2i 



as defined in Eq. (15), showing that light neutrino mixing is no longer unitary. The addi- 
tional possibility of Wl — Wr mixing allows for diagrams with, for instance, two left-handed 
hadronic currents but one left- and one right-handed leptonic current [see Fig. 4(b)], with the 
corresponding suppression factor of tan^ [Eq. (22)]. 

Oz//5/3 processes in the LR model can be categorized in terms of their topology and the helicity 
of the final state electrons; the most relevant diagrams that will be discussed in detail in what 
follows are shown in Figs. 1, 3 and 4 (see Refs. [Ki] for a complete list). Table 1 contains 
a summary of the relevant amplitudes as well as limits on the particle physics parameters 
calculated using the recent KamLAND-Zen limit [ : ] and the matrix elements in Table 2. 
Note that the chiral structure of the matrix element means that the neutrino propagator 
becomes [18] 



^ + TJli 

Pl,r—^ oPLfi 

— m| 



rrii 



— mf 



or Pi 



— 



9 2 ' 



(29) 
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(a) X (b) A 



Figure 1: Feynman diagrams of 01/(3(3 in the left- right symmetric model, mediated by (a) hght 
neutrinos (the standard mechanism Au) and by (b) heavy neutrinos in the presence 
of right-handed currents {A^^J. The diagram with heavy neutrino exchange and 
left-handed currents («4.^^) is the same as diagram (b), with all particles left-handed 
and the replacement V^* -H- S^^. Diagrams with light neutrino exchange and right- 
handed currents are negligible. 



leading to mass or momentum dependence when the leptonic vertices have the same or op- 
posite chirality, respectively, and providing a useful way to categorize the different possible 
mechanisms. In order to give a very rough estimate of the relative magnitudes we denote the 
masses of all particles belonging to the right-handed sector {Mi, Wr and Sr) as i? ^ TeV, 
and those from the left-handed sector as L ~ 10^ GeV (corresponding to the weak scale, or 
the mass of the Wl). The matrices T and S describing left-right mixing can be written as 
L/R, and the gauge boson mixing angle ^ is of order {L/Kf. Note that with this definition 
the order of magnitude of the type I seesaw contribution is my ~ L'^ / R, which is far too large 
in the naive case (without cancellations), but the estimates made above are still reliable. The 
typical scale of momentum transfer is |g| ~ 100 MeV. 

Mass-dependent mechanisms 

In this case the emitted electrons have the same chirality and there are either hght or heavy 
neutrinos exchanged, with mass denoted by rrii and Mj. With both electrons left-handed the 
amplitude is proportional to 

^ + 2 tan e + tan^O E " f) ' (30) 

whereas if both are right-handed it becomes 
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Normal Inverted 




Figure 2: The effective mass (mee) as a function of tlie liglitest neutrino mass in botli tfie 
normal and inverted ordering, with the oscillation parameters varied in their 3a 
ranges [10]. CP conserving (violating) areas are indicated by black lines (blue 
hashes), and prospective values of Yl^i ^^"^ shown. 

Here we have taken into account diagrams with gauge boson mixing at one or both vertices, 
but the most relevant diagrams are: 

• Fig. 1(a), the "standard" diagram, with an amplitude proportional to 

(32) 

where ji?^! — (100 MeV)^ is the typical momentum exchange of the process. The particle 
physics parameter |(mee)| = E^ei"^*! is called the effective mass, and the suitably 
normalized dimensionless parameter that describes lepton number violation is 

M = ^ = ^^<7.1xl0-^ (33) 

nie rrie 

with Uei the (PMNS) mixing matrix of light neutrinos and rrii the light neutrino masses. 
Here and in what follows we give limits on the particle physics parameters 77^; they are 
explicitly defined in Section 3.1.2. The currently allowed [19] regions of the effective 
mass are plotted against the lightest mass in Fig. 2. We will translate this plot into 
half-life in the following section; 

• Fig. 1(b), which is the analogous diagram with purely right-handed currents, mediated 
by right-handed neutrinos. The amplitude is proportional to 

-^L-Oli^Yz'-i^^, (34) 
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(b) As, 



Figure 3: Feynman diagrams of double beta decay in the left-right symmetric model, mediated 
by doubly charged triplets: (a) triplet of SU{2)r and (b) triplet of SU{2)l. 



where (j^Wl) is the mass of the right-handed Wr (left-handed Wl), Mi the mass 

of the heavy neutrinos and V the right-handed analogue of the PMNS matrix U . The 
dimensionless particle physics parameter is 



ei 



< 7.0 X 10" 



(35) 



A diagram not shown in which heavy neutrinos are exchanged with purely left-handed 
currents. The amplitude is proportional to 



A 



C2 t2 
r^2 \^ "^ei ^ ^ 



(36) 



with S c:^ L/R describing the mixing of the heavy neutrinos with left-handed currents. 
The limit is 



C2 



< 7.0 X 10" 



(37) 



Note that the sum in Eq. (36) can be written as 



S] 



(3^ 



which vanishes for negligible Dirac Yukawa couplings. It is also possible to have light 
neutrino exchange with right-handed currents [the term proportional to T in Eq. (31)], 
but this diagram is highly suppressed. 
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Figure 4: Feynman diagrams of double beta decay in the left-right symmetric model with final 
state electrons of different helicity: (a) the A-mechanism and (b) the ?7-mechanism 
due to gauge boson mixing. 



Triplet exchange mechanisms 

• Fig. 3(a) is a diagram with different topology, mediated by the triplet of SU{2)ji. The 
amplitude is given by 

\^WrJ V ^ 



and the dimensionless particle physics parameter is 

< 7.0 X 10- 



= 2^ "4 7^ ^ 7.0 X 10-^ (40) 

Here we have used the fact that the term \/2vRhee is nothing but the ee element of 
the right-handed Majorana neutrino mass matrix M/j diagonalized by V [cf. Eq. (12)], 
with vr the VEV of the triplet 6r and hee the coupling of the triplet with right-handed 
electrons, so that this diagram still indirectly depends on the heavy neutrino mass; 

Fig. 3(b) is a diagram mediated by the triplet of SU (2)^, also present in the usual type II 
seesaw model (without left-right symmetry). The amplitude is given by 

Asl^GI^, (41) 

^L 

which is suppressed with respect to the standard light neutrino exchange by at least a 
factor q^/ml 
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Momentum dependent mechanisms 



In this case the emitted electrons have opposite hehcity, and the amphtude is proportional to 



A 



+ tane + ^ tane + tan^^) E f^-T,*^ - ^.K*^) ; (42) 

the most important diagrams are those involving light neutrinos and two powers of the left- 
right mixing in the pre-factor, i.e., 



The so-called A-diagram in Fig. 4(a), with an amplitude 



~ G 



™WrJ 



1 

UeiT*i - OC 



q R^q 



(43) 



and particle physics parameter 



TUWr 



< 5.7 X 10"^ 



(44) 



Note that this is a long-range diagram with light neutrinos exchanged, with the matrix 
Tg* = 0{M]:)/ Mji) quantifying the mixing of light neutrinos with right-handed currents. 

The ?7-diagram in Fig. 4(b), which also has mixed helicity and light neutrino exchange 
(long-range diagram). This is only possible due to Wl — Wr mixing, described by the 
parameter tan^ [see Eq. (20)]. The amplitude is 

L3 



i 

with particle physics parameter 



OC 



q B?q 



(45) 



\r]ri\ = tan^ 



< 3.0 X 10" 



(46) 



Ref. [jU] gives a detailed explanation of how a complicated cancellation of different 
nuclear physics amplitudes leads to a limit on the r^-diagram that is much stronger 
than the one on the A-diagram. The heavy neutrino contributions to both the A- and 
ry diagrams are further suppressed, being proportional to >S'gjVg*g/M/ [see Eq. (42)]. 
Using the mixing matrices in Eq. (28), the relevant sums become 



s y* 



(47) 



where we have omitted third order terms. This again shows that the left-right mixing is 
a ratio of two scales, M/) and Mr. 
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Using our rough estimates (in terms of L ^ 10^ GeV and R ~ TeV) of the scale of each 
diagram we can now make a naive guess at their expected relative magnitudes. Since the 
mixed A- and //-diagrams in Fig. 4 are of order {L/R)^/q and the purely right-handed short- 
range diagrams in Figs. 1(b) (heavy neutrino exchange and right-handed currents) and 3(a) 
{SU{2)ji triplet exchange and right-handed currents) are of order L^/R^, we expect the mixed 
diagrams to dominate by a factor R^/{Lq) ~ 10^. In the same sense, the amplitudes of the 
mixed diagrams are also larger than the one for heavy neutrino exchange with left-handed 
currents, proportional to I? jR? . The main point is that mixed diagrams should be examined 
more thoroughly, as has been done in the context of inverse neutrinoless double beta decay 
[35] and inverse/extended seesaw [16, 17]. The reliability of the rough approximations made 
here will be discussed in the following section. It should be obvious that the light neutrino 
mass from type I seesaw, my ^ M'^/Mr — L"^ / R cannot be small enough without matrix 
cancellations; the crucial point is that the left-right mixing M£,/Mr ':^LjR can still be large 
in those cases. 

The reliability of these rough approximations can be tested by normalizing the amplitudes 
to the standard contribution, using known bounds on the left-right mixing. The limits on 
the difference of the diagonal elements of the product e^, = [S'S'"'"]q,q, ~ [T'''T]qq, from lepton 
universality [Tl] are 

ee - < 0.0022 , e^, - < 0.0017 , ee - < 0.0039 , (48) 

which give a rather weak bound on the left-right mixing. In the absence of cancellations one 
finds a much stronger bound, namely 

|5..|^|Ti:|^,/^<10-^f^V^', (a = e,/x,r), (49) 



which we apply in the estimates that follow, along with the light neutrino mass scale my 
0.05 eV and momentum exchange |g| ~ 100 MeV. 

For heavy neutrino exchange with right-handed currents [Fig. 1(b)] we have 

^ V Kf TeV \ VTeV\ 



^ ^ ™ \ ^ ^ 8.36 \^\ Hr?^ , (50) 

Ay \mw^) ^ Mimy \mwnJ \ J ' 

whereas for heavy neutrino exchange with left-handed currents [Eq. (36)] the ratio is 

One sees immediately that this process requires cancellations to be enhanced^. However, for 
the A- and 77-diagrams [Fig. 4] we have 

A,^A,^(m^y^ t/.T:.^ <(^y -J^ ^ 2.89 ^ ^ ' ^^^^ 



Ay Ay \mwnJ ^ \^Wr/ y/myMi \t^WrJ \Mr 



(52) 



^This case was also studied in Ref. [31, 33, 34, 52]. 
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Table 1: Summary of relevant mechanisms for Oz//3/3 in the left-right symmetric model, with 
limits on new physics parameters (written in bold face) in each case (see also Ref. [2]). 



mechanism 



amplitude 



current limit 



A-mechanism (Ax) 
?7- mechanism (^r;) 



light neutrino exchange (Au) 
heavy neutrino exchange (Aj^f^) 

heavy neutrino exchange (A^^) G^rriyy^ 
Higgs triplet exchange (Asj^) 



q2 \ ez « 



SI, 



r*2 



S~~ 





UeiT* 







G' 



Gl-l^^^^EiUeiT:, 



0.36 eV 

7.4 X 10-9 GeV"^ 
1.7 X 10-16 GeV-^ 

1.7 X 10-16 GeV-^ 

8.8 X 10-11 GeV-2 
3.0 X 10-9 



where the first inequality comes from the upper limit |^| 



< 



. Depending on the relative 



magnitude of the bi-doublet VEVs Ki and K2, the amplitude Ar, may be further suppressed 
[see Eq. (22)], but this could be compensated for by the fact that M^^" ^ lO^A^^'' (cf. Table 2). 
The main point is that with and Mr around the TeV scale the amplitudes A^^, A\ and 
Ar^ turn out to be quite close in magnitude, whereas the small value for Aj(f^ could still be 
enhanced by cancellations. The mixed helicity diagrams Ax and Ar, can therefore compete 
with the standard light neutrino diagram, even with the stringent limit on T in Eq. (49) 
that connects the left-right mixing to light neutrino mass. Note that in order to arrange 
for cancellations, the Yukawa matrices / and / need to have non-trivial flavour structure so 
that the correct light neutrino mass [see Eq. (12)] can be obtained, since with 0{1) couplings, 
oc Kj, so that Mp, would be near the electroweak scale of 10^ GeV. Assuming that k,2 ^ 1^1 
(see also Ref. [ 1.0]) means that Md — Kif /\/2 and Mi ~ / \/2, so that one has the freedom 
to choose / without affecting the charged leptons. 



3.1.2. Nuclear matrix elements and lifetime 

In order to translate the dimensionless particle physics parameters r^^ into actual lifetimes 
of Oz//3/3 processes for different isotopes one needs the relevant nuclear matrix elements and 
phase space factors. There are various different methods to calculate those quantities and 
most previous studies have focussed on the standard light neutrino exchange mechanism; here 
we attempt to compile a list of the most recently calculated matrix elements relevant to Oz//3/5 
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in the LR model, combining the calculations of various groups. 

We use the QRPA calculation of the matrix elements for the mixed diagrams in Ref. [53] 
(see also Refs. [54, 55]). In their notation, the lifetime of Oz//3/3 can be written as 

+ C4\r]x\^ + + CelvxWVvl cos(V'i - ^2) + Re C2Xn7]x + C^XrT]^ } , (53) 

where the coefficients Ci are combinations of matrix elements and integrated kinematical 
factors, Gqi is the usual phase space factor and ipi are complex phases. The parameters 
Xl {Xfi) include all processes in which the final state electrons are both left-handed (right- 
handed), i.e. 

Xl = Mtvu + -M'^'r^^^ + M"j:;vs, , and Xr = Al'^'r/^^ + A^'^^^r/^^ , (54) 

with T]Sj^ the LNV parameter associated with Eq. (41). In Eq. (53) we have omitted the 
interference term XlXr, which is suppressed due to the different electron helicities (e^e^ 
vs e]^e^); interference terms with final states in which at least one of the electrons has the 
same helicity have been included. The matrix elements Ai'^"^ and A^'^ include Fermi and 
Gamow- Teller contributions. 

Ref. [50] presents an improved calculation of the phase space factor Gq'( for the light neutrino 
exchange mechanism, taking into account the finite nuclear size of the Dirac wave function as 
well as electron screening effects and angular correlations. The factor is slightly lower, with 
the difference becoming more marked for heavier nuclei. The coefficients Ci {i = 2,3,4,5,6) 
depend on different phase space factors [53, 57]; here we assume those factors are reduced 
by the same percentage as G^'^. More recent calculations [58-60] of light and heavy neutrino 
matrix elements include the Gamow- Teller factor A^gt the relevant matrix elements Ai^" 
and Ai^. For consistency of notation, we make the following definitions 

I I (55) 

which allow us to write the lifetime in Eq. (53) as 

+ |Al°''n?7Ap + lAfJ'^n^^r,!^} + interference terms. (56) 

The corresponding matrix elements are reported in Table 2 and will be used in the analysis 
that follows. The range of values comes from the fact that different calculations have been 
used. Note that we have used the new phase space numbers to calculate limits. 

In the limit of type II seesaw dominance, the expression in Eq. (56) will simplify considerably, 
whereas with type I seesaw dominance all six terms should be considered (we neglect the 
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Table 2: The phase-space factor Gq^ [53, 56] and the matrix elements for light {Ai'l'^) [61] and 
heavy {Ai'j!f) [59, 60] neutrino exchange, and for the A- and 77-diagrams [53, 55], for 
different isotopes, for = 1-25 and tq = 1.1 fm, corresponding to Eq. (56). 



Isotope 


[10-1^ yrs-i] 
(old[ ]) 


[10-14 yrs-i] 
(new [ ]) 


Ml" 




Ml'' 




^^Ge 


0.793 


0.686 


2.58-6.64 


233-412 


1.75-3.76 


235-637 


82Se 


3.53 


2.95 


2.42-5.92 


226-408 


2.54-3.69 


209-234 




5.54 


4.13 


2.43-5.04 


234-384 


2.85-3.67 


414-540 


i36Xe 


5.91 


4.24 


1.57-3.85 


160-172 


1.96-2.49 


370-419 



contribution stemming from the left-handed triplet Sl, which is suppressed by light neutrino 
mass and ^5-- = C(TeV)). We use the notation [T-jy'g]^ {k = u, N^\6ji, X,ri) to refer 
to the lifetime stemming from one particular diagram. Figure 5 illustrates the variation of 
the lifetime [T'{/'2]i' with lightest neutrino mass, for the Oz//3/3 of ''^Ge and using both the 
smallest and largest matrix element {M^^ = 2.58); comparison with Fig. 2 shows that the 
lifetime is obviously just the inverse of the effective mass, with various numerical pre-factors. 
The variation in Ml'^ can bring the minimum allowed lifetime down by roughly one order of 
magnitude, and affects the translation of the KamLAND-Zen limit into a limit on the half-life 
for ''^Ge, as discussed in the figure caption. 



3.2. Charged lepton flavor violation and dipole moments 

Although small active neutrino masses "GIM suppress" charged lepton flavor violating pro- 
cesses by a factor of (Am^/m^^)^ ^ 10^^° {Am\ is the atmospheric mass squared difference), 
the existence of heavy right-handed neutrinos and Higgs scalars allow the LFV decays ^ ^ e'j 
and /i — )■ 3e as well as — )■ e conversion in nuclei to occur at rates observable in current ex- 
periments. Those decay rates will be proportional to similar combinations of mass and mixing 
parameters as the Oi//3/3 amplitudes, thus providing complementary constraints. Defining 

= r(/i- ^ e-z/^Pe) and Feapt = r(/i- + A(Z, N) ^ + A(Z - 1, N + 1)), (57) 

the relevant branching ratios 

_ r{f,+ ^ e+7) 

_ rif^- + AiN,Z)^e' + AiN,Z)) 

^ capt 

_ F(/i+ ^ e+e-e+) 

DiXfj,^3e = ^ , 
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Figure 5: The standard light neutrino contribution to the Oz//3/3 half-hfe of ^^Ge plotted against 
the lightest light neutrino mass, using 3a ranges of the oscillation data from Ref. [49]. 
Shaded regions (dotted lines) are for the smallest (largest) NMEs from Table 2. The 
grey shaded region is excluded by the KamLAND-Zen experiment, the horizon- 
tal dashed (dashed-dotted) lines show the planned sensitivities of the GERDA [62] 
experiment, with 40 kg (1 ton) of isotope. The Heidelberg-Moscow limit [63] is in- 
dicated by a horizontal (red) dotted line. The dotted black line shows the smallest 
value of the KamLAND-Zen limit, with the NMEs for ^^Ge and ^^^Xe varied in their 
allowed ranges from Table 2. 



are constrained at 90% C.L. to 

BR^^e7 < 5.7 X 10"^^ [ ], R;^^, < 7.0 X lO'^-'^ [ ] and BR^_,3e < 1.0 x lO"^^ [ ] 
by experiment. 

The amplitudes for LEV decays in the LRSM receive contributions from (i) right-handed 
gauge bosons and Higgs triplets, suppressed by {mwj^/mwj^)^; {ii) left-handed gauge bosons, 
suppressed by ~ \Mj^M^^\'^ and {Hi) processes with Wl — Wr mixing, suppressed by ^M^M^^ . 
Terms proportional to are expected to be small and are neglected here. All of the possible 
channels are in some way related to the right-handed neutrino mass, either directly as a virtual 
particle in the loop or indirectly since the couplings of the triplets to leptons are proportional 

to Mr^. 

A detailed calculation of the LEV decay widths and branching ratios in the LRSM has been 
performed in Ref. [f)7], where the results have been obtained by expanding to leading order 

''The assumption of a discrete left-right symmetry means that the exchange of left-handed triplets is also 
related to right-handed neutrino mass, see Eq. (A-5). 
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in the ratios M^i/Mr and k^/vr, and thus ignore any effects of left-right mixing. The results 
are (see also Refs. [68, 69]) 



o T-) triplet 1 
^%^3e = o 



h^eKe 



2 / m^^ m^^ 



4 4 / ' 

<5++ "^5+^ 
"l "r 



(59) 



for the tree-level process — )■ 3e and 



1 T V \ ^ 

BR^^e7 - 1-5 X 10-^ ( — , (60) 



1 TeV \ / , ^s++ 



2 



-8xlO-|W^^— j "[log-|^| , (61) 

for the loop-suppressed decays fi ^ ej and /x — )■ e conversion, where the expressions are 
simplified by assuming the "commensurate mass spectrum" Mj ~ — fTT's++ — f^s++ — 

'm^+. The parameters h and g\i^ are defined to leading order in the ratio Mo /Mr by 

R 



/.., = X,^..V^,.^ = ^ and ^ifv = l^^,.K.^;;^ 



-RJ lie 



(62) 



assuming manifest left-right symmetry (i.e., a discrete parity symmetry, see the Appendix). 
If one assumes that logarithmic terms [see Eq. (A-14)] from doubly charged Higgs diagrams 
dominate and that no cancellations occur amongst the LFV parameters {\guv\ — \h^eKe\)^ '^^^ 
expects BR^_^3e to be roughly two orders of magnitude larger than R^_j.e for (9 (TeV) Higgs 
triplet masses [u.]. Thus in this simplified case the limits on /x — )■ 3e will confine the model 
parameter space the most. 

However, with right-handed neutrinos around the TeV scale the left-right mixing could 
be enhanced, so that the usual type I seesaw contribution to LFV processes should also be 
considered. Those have been calculated in Refs. [33, 70-77]. Since the LRSM is effectively a 
type I-I-II seesaw model one needs to take into account LFV processes mediated by both heavy 
neutrinos and Higgs triplets, effectively allowing for interference between different amplitudes. 
Ref. [77] has presented the full expressions for ^ — )■ 67; we include type I seesaw terms in 
the expressions for /i — )■ 3e and yU — )■ e conversion, in the former case including possible 
interference between loop and tree level diagrams. Detailed expressions for the decay widths 
including form factors and loop functions can be found in the appendix; we summarize the 
most constraining processes here. In our parameter scans in the type I dominance case we 
take into account all relevant contributions. 

It turns out that the most important constraints on the mixing S ~ M^/Mr come from 
/i — )■ 67 and /x — e conversion. In both cases the constraint is roughly 

s;,s,,7{x;) ^ s;,s^ < lo-^ , (63) 
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where we take the loop function T{xi) to be of order one. This approximation is not always 
valid for very large right-handed neutrino masses, in which case J^(xj) ~ ln(M?/m^^), but 
since the mixing scales with 1/Mj the rate will vanish in the decoupling limit [ ]. The 
loop-suppressed decay rate of yU — > 3e (with heavy neutrinos exchanged) depends on the 
same parameters as /i — t- e conversion, but the limits are weaker in this case: the bound 
BR^^Se < 1.0 X 10~^^ can be roughly translated into 5'*jS'gj ^ 10^^. These constraints come 
from diagrams with left-handed currents and left-right mixing, i.e. the terms proportional 
to in Eqs. (A-9), (A-14), (A-15) and (A-17), so that there is no other dependence on the 
heavy particle masses besides from the loop functions. Another interesting constraint comes 
from — )■ 67 diagrams in which gauge bosons mix: the chirality flip occurs within the loop, 
leading to a direct dependence on the Dirac mass matrix instead of the muon mass [Eq. (A-9)] , 
in a similar way to the mixed diagrams in Oi//3/3 (see also Refs. [74, 77, 78]). This enhances 
the contribution of mixed diagrams to yU — )■ 67 by a factor SMp/^/m^ ~ Mo/frif^, so that the 
product of the mixing angle C, and the fie element of the Dirac mass matrix is constrained to 
be 

l^^u(l^) ^0.2GeV. (64) 

In addition, the experimental limit of \de\ < 10^^^ e cm [ ] on the electric dipole moment of 
the electron [see Eq. (A- 10)] constrains the ee element to be roughly 



lm{[Mn]eee'''} (^-^^ ~ 0-02 GeV, 



(65) 



which also depends on the phase a. These limits effectively constrain the ?7-diagram in 
Fig. 4(b). 

One might also expect large left-right mixing to allow loop-suppressed (type I) contributions 
to /i — 3e to compete with the tree level triplet (type II) contribution. The full expression 
is given in Eq. (A- 12), and the condition for comparable magnitudes of type I and type II 
contributions is roughly 

'5TeVV flMM* 



assuming m^++ = nT's++ = tt^s++- Thus for TeV-scale Wr the bound on S"^ in Eq. (63) 
means that one needs right-handed neutrinos around the electroweak scale for the type I loop 



contribution to be competitive in yU — )■ 3e decay. 



3.3. Collider physics 

Before concentrating on the Oi//3/3 amplitudes we briefly discuss the role of collider physics in 
studying the LRSM. Colhder searches provide a complementary probe of the parameter space 
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Figure 6: Comparison of the limits in M^v^ — rriw^ parameter space from CMS and from the 
KamLAND-Zen hmit on Oz//3/3. The hmit of 1.9 x lO^^ yrs on the Oz//3/3 half-hfe 
of ^^^Xe means that all points to the left of the solid black line (dashed red line) 
are excluded, for iKip = 1 = 2/3), where we assume only heavy neutrinos 

contribute to 0^(3/3, i.e. only [Ty2]j^ifi)- The shaded region is excluded by CMS at 
95% C.L. [80]. 



of the LRSM: the right-handed W boson and right-handed neutrinos can be produced in pp 
collisions at the LHC via [81] 

PP^Wr + X ^ N, + i + X , (£ = e,/i), (67) 

followed by the decay into like-sign dileptons and two jets, i.e. 

Wr iiNe ^ £i£2l^5 ^ iii2qq' ^ iihjj , (68) 

which for the i = e case is equivalent to the Oz//3/3 diagram in Fig. 1(b). The CMS collabora- 
tion looked for this signature in both 7 TeV [10] and 8 TeV [80] data, where the integrated 
luminosity was 5.0 fb^^ and 3.6 fb^^, respectively. Their analysis was simplified by assuming 
negligible mixing (^ ^ 0) between gauge bosons and between heavy neutrino mass eigenstates 
(V ~ l), so that the final states are either both electrons or both muons. ATLAS studied the 
same process with 2.1 fb~^ of data from 7 TeV collisions [ ], and in addition examined the 
case of maximal mixing between the first two heavy neutrino mass eigenstates. 

As a simple illustration of the complementarity of the different data sets we plot the limits 
from the latest CMS data as well as from the KamLAND-Zen OuPP experiment [47] in the 
M^r^ — mwj^ parameter space in Fig. 6, using two different values for the mixing Vei- Here 
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one assumes that only one heavy neutrino flavour N^, ^ A^^i is accessible, so that the LNV 
parameter in Eq. (35) simply becomes |?7^^| = mp(mvKi, /""^Wj;)^ I KiP/^i- 

It is also possible to probe the couplings hap of Higgs triplets to leptons [see Eq. (5)] with 
collider searches. The latest results from ATLAS [45] give the exclusion limits Tn^±± > 409 GeV 
and m^±± > 322 GeV for e^e^ final states and assuming a branching ratio of 100% to each final 
state. In order to compare these results to the Oz//3/3 bounds one needs to take into account 
the other decay modes of doubly-charged Higgs scalars into gauge bosons and singly-charged 
scalars. An analysis in this direction was performed in Ref. [82], and the results depend largely 
on the mass spectrum of the different components of the Higgs triplets /j. 

4. Oz//3/3 amplitudes in the seesaw limits 

In the most general case the light neutrino mass matrix 

m, = ML- MdM^^MI , (69) 

receives contributions from [see Eq. (13)] both the left-handed triplet (type II seesaw) and 
the heavy right-handed neutrinos (type I seesaw), making quantitative studies of the Oi//3/3 
amplitudes difficult. Here we focus on the two extreme cases of type II and type I dominance; a 
complete study is beyond the scope of this work. In the former case one sets the Dirac Yukawa 
couplings to zero, in the latter one assumes that the triplet VEV vanishes, i.e., vl = 0. The 
simpler case of type II seesaw dominance is dealt with first. 



4.1. Type II seesaw dominance 

With the approximations mentioned above, the lifetime in the limit of type II dominance is 

[^i%];pe n = G'oi { |-M°f + l-M^f \vk + ^^.1'} ; (70) 

by neglecting all Dirac Yukawa couplings we drop all terms proportional to Md, i.e., those 
with left-right mixing. We are left with only heavy neutrino [Fig. 1(b)] and triplet exchange 
[Fig. 3(a)] in addition to the standard diagram [Fig. 1(a)] (the amplitude Aj^^ also vanishes, 
being proportional to M^)- As discussed above, the interference term is suppressed, since the 
final state electrons in Fig. 1(a) are left-handed whereas those in Fig. 1(b) are right-handed. 

In the case of type II dominance, the right-handed neutrino mass matrix can be expanded 
as [36] 

j\^type II ^ Vr ^ n'ihDTn'^h^ - k'\_— — {h om'^'^ h D)m~^ {h om'^ h of + . . . , (71) 
and since we neglect Yukawa couplings (/id ~ 0), 

= ' (72) 
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Figure 7: Same as Fig. 5, with the grey shaded areas forbidden by the LFV constraint 
BR^_j.3e ^ 10^^^, for different values of the Higgs triplet mass, with = 3.5 TeV 
and the heaviest right-handed neutrino Mheavy = 500 GeV. 



which simplifies the analysis considerably: the light and heavy neutrino spectra are propor- 
tional to each other, and V = U, up to an overall complex phase. In addition, both U and 
V become unitary in the limit that = [cf. Eq. (28)]. These assumptions were first used 
in Ref. [11] to quantify the heavy neutrino contribution to neutrinoless double beta decay; 
the triplet contribution to Oz//3/3 was neglected since the constraint from /i — )■ 3e leads to 
Mn/msj^ ^ 1 over a large part of parameter space. We recast this idea in terms of actual 
lifetimes, and show explicit regions in parameter space where the limit from BR^^3e comes 
into play. Replacing V with U in Eq. (35), the dimensionless LNV parameter corresponding 
to heavy neutrino exchange with right-handed currents (oc [M^^jee) can now be written as 

[WNO-m,^---J |f/,i| +— e +|f/e3|e j—, (73) 

[^ivJio-m,^---J |f/,i| +|f/,2| e +— If/esle J—, (74) 

for normal and inverted ordering, respectively, where a and (3 are Majorana phases. Similarly, 
the branching ratio for /i — )■ 3e in Eq. (59) depends on the product of the ee and fie elements 
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of h = Mji/mwn, with 

[Mr]'! = i—U.iU,i + f/.2f/p2e^" + —U^.U^^eA . (75) 



We assume ^s++ — ^s++ what follows. 

Following Ref. [ ], by fixing = 3.5 TeV and the heaviest right-handed neutrino mass 
Mieaviest = 500 GeV, the three contributions can be plotted against the hghtest light neutrino 
mass (see Figs. 7 and 8). It is clear that the right-handed contribution [Ti/2]~(fl) [Fig. 8(a)] 

is proportional to the inverse of Mr, whereas the triplet contribution [T'i/2]5jj^ [Fig. 8(b)] is 
proportional to Mr, and looks similar to the standard lifetime [Ti/2],7^ (Fig. 7), since m^, oc Mr 
in the type II limit. For [T'i/2]"!^) , the inverted ordering can have infinite lifetime (zero effective 

It 

mass), whereas the normal ordering cannot, so that the roles are reversed with respective to 
the standard case. In each plot we indicate the regions excluded by the limit on /i — )■ 3e 
for different values of m^++: in the normal hierarchy the constraint only comes into play 
when the lightest mass is larger than about 0.01 eV, whereas in the inverted hierarchy the 
whole parameter space is affected"*. In the case of the light neutrino and triplet contributions, 
the only areas still allowed correspond to the largest possible value of (rriee), i-e., when both 
Majorana phases are close to zero. 

Figure 9 shows the total half-life, with all three contributions included. The chosen value of 
771^++ affects not only the LFV constraint but also the resulting half-life, due to the dependence 
of the triplet contribution on this quantity [Eq. (39)] . The black dotted lines show the half- life 
without the triplet contribution, and it is evident that the addition of the triplet part can 
shorten the half-life by several orders of magnitude, bringing it within reach of the GERDA 
experiment. There also exist regions where the lifetime can be longer, due to cancellations 
between the t]^^ and r/^^ contributions. The key point here is that the triplet contribution can 
still be allowed for certain values of the Majorana phases, even with the LFV constraint, thus 
enhancing the total amplitude for OuPP. This enhancement obviously depends on the triplet 
mass, so that if fns++ ^ 5 TeV we recover the results of Ref. [11]. 

4.2. Type I seesaw dominance 

In the limit of type I seesaw dominance all the terms in Eq. (56) must be considered (we 
neglect the small contribution from 775^, as discussed above). This leaves us with six con- 
tributing diagrams: (i) "standard" light neutrino exchange {r]^)', (ii) heavy neutrino exchange 
with left-handed currents {Vnr)'^ (^^^) heavy neutrino exchange with right-handed currents 



^Our results agree with Fig. 2 of Ref. [11], which shows that Mhcavy/'7i|5++ ^ 0.1 in the inverted ordering for 
aU light neutrino masses, which in our case would correspond to 1^^++ ~ 5 TeV. 
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Figure 8: The contribution to the OuPP half-hfe of ^^Ge from (a) heavy right-handed neutrinos 
and (b), right-handed Higgs triplets, plotted against the lightest light neutrino mass, 
with = 3.5 TeV and Mheavy = 500 GeV. In plot (a) the grey shaded regions 

are excluded by LFV constraints, for different values of rng++, in plot (b) ^s++ ~ 
tt^Wr = 3.5 TeV. Experimental limits are explained in the caption of Fig. 5. 
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Figure 9: The total Oz//3/3 half-life of ^^Ge including light neutrino, heavy neutrino and triplet 
contributions, plotted against the lightest light neutrino mass, with mv^^ = 3.5 TeV 
and Mhcavy = 500 GeV. The solid, dashed and dashed-dotted lines show the allowed 
regions that satisfy BR^^2,e < 10^^^ for equal to 1, 2 and 3.5 TeV respectively; 
the black dotted lines enclose the regions allowed if one neglects the triplet contri- 
bution and LFV constraints. Experimental limits are explained in the caption of 
Fig. 5. 

{rj§^); {iv) light neutrino exchange via the A-diagram (//a); (v) light neutrino exchange via 
the 77-diagram {rjri) and {vi) right-handed triplet exchange {r]Sj^). There are also interference 
terms [see Eq. (53)], and distinguishing the different contribution becomes difficult. The con- 
tributions (iv) and (v) are usually neglected in the literature, but can actually be significant, 
as we have shown in the rough estimates above. 

4.2.1. Parameterizing the relative magnitudes 

In order to quantify the six contributions one needs more information about the right-handed 
sector, specifically the right-handed mixing matrix Vr and the mass spectrum Mj (i = 1, 2, 3) 
of right-handed neutrinos. The right-handed mass matrix M/j appears in the amplitudes A'^^, 
Aff^, Asjf, Ax and Arj, and in the case of type I seesaw dominance can be expanded as 

The leading term is a matrix product containing the unknown Dirac mass matrix, so that the 
simple relations in Eq. (72) no longer hold and one needs a different approach. The authors 
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of Ref. [ I i] simplify the analysis by assuming that (i) the Dirac mass matrix is diagonalized 
by Vr and (ii) the three Dirac Yukawas are equal. This scenario is very restrictive; another 
approach would be to insert an ansatz for the matrix of Dirac Yukawa couplings ho. Often one 
uses the condition M„ ~ = n+hf), which holds at the GUT scale in 5*0(10) models [83]. 

More generally, the Dirac mass matrix can be parameterized using the so called top-down 
or "VL-parameterization" 

Md = uIMdUr, (77) 

where Ul and Ur are arbitrary unitary matrices and = k+ diag(/;,i, h2, h^). In the LRSM 
type I case, M^) has 18 parameters and Mr has 12 parameters, so that the left-right mixing 
MjjM^^ depends on 30 parameters, making it difficult to learn anything from a parameter 
scan. If we assume a discrete parity (charged conjugation) symmetry, then becomes 
hermitian (symmetric) thus reducing the number of parameters by 6. However, it is still 
numerically difficult to find Dirac mass matrix structures that give large enough left-right 
mixing. One way is to start from a specific matrix structure in M^) that gives zero neutrino 
masses, and introduce small perturbations (see Refs. [29, 84]). 

It is indeed possible to scan the entire allowed parameter space using the orthogonal pa- 
rameterization [x~i], where the Dirac mass matrix is written as^ M^i = iV,yml,^'^OMj^'^V^ , 
with 00^ = 0^0 = 1 and the diagonal matrices fhi, = diag(mi, m2, ms) and Mr = 
diag(Mi, M2, M3). However, due to the large number of unknown parameters (6 in O, 3 
in Mr and 9 in Vr) it is difficult to learn anything. In addition, the 0-matrix approach does 
not allow one to define a symmetric or hermitian Dirac mass matrix in a simple way. 

An alternative method is to go to the basis where Md is "diagonal", so that the light 
neutrino mass matrix is given by 

ml = -MoM'^-^Md , (78) 

with M'^~^ = UrM^^UI. In essence one has rotated the left-handed neutrino fields hy Ul 
[cf. Eq. (77)]. After diagonalizing m'^, by the unitary matrix X^, i.e. m'^, = XifhvXj^, the 
neutrino mass matrix in the fiavour basis is 

= -V^Xl [MnM'jf^Mn) XlVj = -Ul (^MdM'j,-'Md) Ul , (79) 

where is the light neutrino mixing matrix [Eq. (16)] defined by V^, = Ul^X^. Numerically, 
this means one needs only fit the mass eigenvalues after diagonahzing Eq. (78), decoupling 
the PMNS mixing parameters^. The authors of Ref. [52] used this approach to find matrix 
structures that could enhance the amplitude for double beta decay mediated by heavy sterile 
neutrinos [Aj^^), albeit without right-handed currents. In our case those same structures will 

^Note that in the left-right model we cannot rotate to a basis where Mj^ is diagonal without affecting the 

right-handed charged current. 
''This approach is discussed in Ref. [86]. 



26 



also enhance the amphtudes for the A- and r^-diagrams and influence the LFV branching ratios. 
However, one cannot recover the non-trivial mixing Vr in the right-handed sector simply by 
diagonalizing Mjj"\ Defining M;^^^ = X^M^^Xjj means that 



Vr = uix 



^0) 



so that the only way to find Vr is to invoke the symmetry (hermiticity) of Md, which gives 
Ur = Ul {Ur = Ul)- The right-handed mixing is then 

Vr = uIXr = V,xIXr or Vr = UIXr = V*XTX 

whereas the left-right mixing (in the flavour basis) is 



L^R ) 



MnM^i = U\Mr>M''^Ur or MnM^^ = U\Mr>M''^m 



12) 



for symmetric or hermitian M^, respectively. The expression [cf. Eq.(38)] characterizing the 
diagram with heavy neutrinos and left-handed currents is 



L ■ 



^3) 



The corrected forms of U and V used for calculating Oz//3/3 amplitudes and LFV branching 
ratios can be found from Eq. (28), but in our case the terms second order in R M^M^^ 
make little difference. 

It has also been shown "] that if the Dirac mass matrix is symmetric, there are only 



2^ = 8 discrete solutions to the seesaw equation, given by Md = i\J MyM]^^ Mr. We have 
checked that it is possible to obtain large left-right mixing solutions that are consistent with 
this formalism. In that case half of the eight solutions give large mixing, whereas the other 
half give small mixing. 

4.2.2. Numerical example 

In the most general case, one should solve the condition M^M^^ M]^ = 0, and it turns out 
that in the basis in Eq. (78) this equates to [52] 

/O l\ 



Md oc diag(0, 0, 1) and M^j oc 



4) 



1 1 

VI 1 ly 

Inserting small parameters instead of zeros leads to non-zero light neutrino masses, with the 
spectrum depending on any hierarchies introduced in Md and Mr. One particular example 
(from Ref. [O-]) is 

fd e f\ 



Md = /t+diag(ae^ be, c), 



R 



~ M-^ 



g he 



15) 
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Figure 10: Contribution to the Oi//3/3 half-life of ™Ge from the A- and //-diagrams plotted 
against the lightest light neutrino mass, for symmetric M/j. The standard con- 
tribution is indicated by the region outlined in black, and the dashed and dotted 
horizontal line correspond to the limits from Eqs. (44) and Eq. (46). 



which leads to nonzero lightest neutrino mass. With all coefficients a, 6, c etc. of order one 
one needs |e| = (9(10^^) in order to get the correct mass for active neutrinos with the matrix 
textures in Eq. (85). Inverting M^~^ would give a matrix with small (1, 1), (1,2) and (2, 1) 
entries, but since M'j^ ^ = [UrM^^UI), the matrix Mr can have large entries everywhere, 
which can enhance the LEV amplitudes. This is simply a manifestation of the fact that one 
cannot go to a basis where the right-handed neutrinos are diagonal without affecting the right- 
handed current, which is different to the conventional case. Eor our parameter scans we set 
fTiWn = 3.5 TeV and 'm^++ = 5 TeV and vary the gauge boson mixing angle in the range 
10^^ < ^ < 10~^ , otherwise it would be difficult to evade the constraints from /i — )■ 67. The 
magnitudes of the complex parameters a, b, c etc. are varied in the range [0.1, 1.0], and |e| in 
the range [10~^^, 10~^]. The phases are taken to be between and 27r, and k+ = 174 GeV and 
M = 1 TeV are fixed. Erom Eqs. (12) and (21) the relation Mr = ^m\Y^h ~ Sniyy^h holds, 
which we used to check perturbativity of the coupling h. 

One expects the different half-life contributions to have similar orders of magnitude, since we 
are exploring the fine-tuned region, so that the amplitudes A^^, Ax and Ar,, which all depend 
on the left-right mixing, are enhanced. We plot the half-lives for the amplitudes A\ and Ar^ 
in Eig. 10 and the half-lives corresponding to heavy neutrino exchange, i.e. the amplitudes 
A^ , A% and Asji in Eig. 11, in both cases for a symmetric Dirac mass matrix. In each case 
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Figure 11: Contribution to the Oz//3/3 half-life of ^^Ge from heavy right-handed neutrinos, with 
left- and right-handed currents (A^^), for symmetric Mjj. The standard contri- 
bution is indicated by the region outlined in black, and the dashed horizontal line 
corresponds to the limit from Eq. (37). 



the usual light neutrino contribution is shown for comparison, and one can see that there are 
regions of parameter space in which the A and t] contributions dominate over the light neutrino 
contribution. Remarkably the r] contribution can still be sizable, even with such small values 
of ^: this is largely due to the larger value of the matrix element Ai^^" (cf. Table 2). The 
lightest mass could be smaller if the parameters a, b, c were allowed to be smaller than 0.1, 
although in the normal ordering case the LFV constraints in general favour larger values of 
'flight- This is largely due to the fact that lighter mught means larger Mr, which enhances the 
decay rates for /i — )■ 3e and other LFV processes. In addition, it turns out that b and c need 
to be small in order to keep the left-right mixing small enough, since the rotation matrices in 
Eq. (82) can lead to large entries in the (1, 1), (1,2) and (2, 1) positions of M^M^\ which 
enhance LFV processes. 

In order to ascertain whether one diagram might dominate over another it is interesting 
to look at the ratios of different half-lives, which has the added advantage that uncertainties 
in NMEs will drop out. In Fig. 12 we show the ratios of various half-lives to the standard 
half-life, calculated for the example texture. Here it is obvious that the A-contribution can be 
larger than the light neutrino contribution. 
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Figure 12: Ratio of half-life contributions, [T^j2]\l[Ty2]y and [T^j2]Nj^l[T^j2]^. 

5. Conclusion 

In this paper we have investigated the interplay of neutrinoless double beta decay and charged 
lepton flavour violation in the context of the left-right symmetric model, paying particular 
attention to those Oz//3/3 diagrams usually neglected in the literature. In the case of pure type II 
seesaw we have shown that the triplet contribution to Oz//3/5 should not be neglected for all 
light neutrino masses. For pure type I seesaw there exist regions of parameter space in which 
all diagrams can have similar orders of magnitude, which makes distinguishing the leading 
contribution difficult. In particular, the momentum-dependent A-diagram can be larger than 
expected. As we have shown, the bounds from lepton flavour violating decays complement the 
study of lepton number violation, and can be used to further restrict the parameter space. A 
comprehensive study should include the type I-I-II case, which we leave for future work. 
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Appendix 



A. Details of lepton flavour violating expressions 

Here we give details of the different contributions to lepton flavour violating processes. 



A.l. Lagrangians 8i couplings 

LFV decays proceed via the charged current in Eq. (23), which we repeat here for convenience, 
as well as the couplings of the charged components of Higgs triplets to lepton doublets in 
Eq. (5); the relevant terms are (with hi = = h) 

9 
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+ h.c. 



(A-1) 



Rotating the fields to the physical basis gives 

,lep 9 



C'^l = [hl^KLULiWr, + ^e^'-W,-^) + enYKRnl{-^e-''^Wr, + W,-^)] + h.c. 



c 



Hi 



Cc±± 



9 
9 



Til 



(A-2) 



where we have used Eqs. (12), (17), (20) and (21), with 

n-L,R = {VL,R) ^R— ^R ^L,R-[^LRI 



R! ~ "L^R^ 

rnwa 



(A-3) 



and My = diag(mi, 1712, m^, Mi, M2, M3). The LFV parameter is 

2 



L,R 



hi oh 



l,r"'L,r 



Mu 



^R ^L,R 
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In the manifest left-right symmetry case (discrete parity symmetry), V[ = V^, so that these 
expressions become [67] 



h = hL = hR = K^^^Kj^ , and gu^ = g^^^ = g^^ 
mwji 



K 



R 



(A-5) 



In our case we take the charged lepton mixing matrices to be diagonal so that all processes 
depend on a combination of the mixing matrices S and V [see Eq. (62)], depending on the 
helicity of the different particles. 
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A. 2. Decay widths and branching ratios 

The effective Lagrangian for /i to e conversion can be written as 



4(47r)2m2 



a 



w 



J2 {ei, [WIPl + W^Pr] /i q^q] + li.c, 



(A-6) 



with cr^i^ = |[7m' '^A the form factors and W^'r- The full matrix element for /x — )• 67 
is given by 



iM{fi — )■ 67) 



eavi/ 



(A-7) 



with the anapole and dipole form factors F2 r and ^ defined in Eqs. (A-14) and (A-9). 

The on-shell decay /i — )■ 67 only receives contributions from the ^ terms, the branching 
ratio turns out to be 



where 
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with Xi = {Mi/mw^Y ^ yi = {Mi/mw^Y and the loop functions G\2{.x) defined in Eq. (A-26). 
In addition, the electric dipole moment of charged lepton ^ct (a = e, r) is given by [15, 77, 87] 

" 3 

■Im 5^5,,\4iee*°G^(a;,)M, , (A-10) 
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which is similar to the mixed diagram contribution in /i — )■ ej. 

The tree level contribution to — )■ 3e in Eq. (59) can be rewritten as 
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to be compared with the loop-suppressed type I seesaw contribution given by [72, 88] 
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The interference terms between triplet exchange and gauge boson mediated loop and box 
diagrams are 



BR 



triplet+type I 
'/i->3e 
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Re [2sl,T*F2 + ^s'^T^Gl + T^B^T + T*FI^{1 - 2^)] (a_13) 



Re [2sl,T*Fl + Asl,T*Gl + T^B'^jl^' - 2s^T*F^^] \ , 



where T = h^eh*^,^, and hap is defined in Eq. (62). Note that the triplet term effectively has 
the same structure as the box contribution (after Fierz transformations, see Ref. [89]), so we 
expect it to interfere with the other amplitudes in the same way. 
The form factors for off-shell photon exchange are 



Fl = Y.\sUS^^F,{x,)-V,,V:,y. 
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where the logarithmic term is a simplified version of the usual triplet loop function [ ], since 
we take the doubly charged scalar mass to be much larger than the charged lepton masses 
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{''^Sl R ^ ^e,fi,T)- The Zi-boson exchange terms^ can be expressed as 

3 

+ (S^S*),, [Gz{x„x,) - Gz{0,x,) - Gz{0,x,)] + {S^ S),,Hz{x„x,)} , 

3 n _ 9 2 2 (A-15) 

- E ^-^el -^^^ (FziVi) + 2Gz{0, y.) - f 

2 2 

H — ^DziVi, Xi) H ^DziVi, Zi 

where Zi = [Mi/mu^Y; the box diagram form factors are^ 

3 

B'lT = -2 E [^Xbox(0, X,) - Fxbox(0, 0)]} 

i=l 

3 

+ E ^U^ej [Fxbox(2;i, Xj) - Fxbox(0, Xj) - Fxbox(0, Xi) + Fxbox(0, 0)] 

~^ 'S'gj5'g^(jrbox(Xi; Xj, ) 

(A-16) 

^^^r = -2^^ E {^M^^e* [i^Xbox(0, yd - i^Xbox(0, 0)]} 



^l^fl i=l 

3 

+ E ^.^^e} {-2K,K: [^Xbox(l/.,%) - -Fxbox(0, - -Fxbox(0, Ui) + -Fxbox(0, 0)] 
+ K*KjG'box(l/i, %■,!)} , 



(A-17) 



for purely left- and right-handed contributions and 

.2 3 / o 



^Lr" — ^ E '^A'i'^eiKjKj^box ( Xi,Xj, 2 ^ ) ' (A-18) 

^ "^V^fl i j=i V ^Wr / 

^^r = E ^..K*^:.^e,Gbox (x., X,, ^) , (A-19) 

^Ri,j=l ^ ^WrJ 

for diagrams with mixed helicity. The loop-suppressed amplitudes with right-handed currents 
contain the 0{1) mixing matrix V as well as the additional suppression factor of (mw^/mvi/^)^; 
without the enhancement from large left-right mixing (in S*), we expect those contributions to 
be much smaller than the tree level one in Eq. (59). The mixed left-right box contributions 



We ignore terms from the exchange of the heavier Z2 boson. 
^We neglect terms proportional to 
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come from an effective four fermion operator, as is the case in kaon mixing [40, 91, 92], with 
a factor of 1/2 coming from the Fierz transformation of a scalar to vector contribution (see 
Ref. [88]). 

— e conversion in nuclei is similar to — j- 3e and receives contributions from the same 
loop and box diagrams. The — e conversion rate is given by [67, 75, 76, 88] 
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are composite form factors. Note that the expression in Eq. (A-20) is derived by approximating 
all interactions to be point-like and taking the proton and neutron densities to be equal. In 
this case the wavefunction overlap integrals D and ^(P'") calculated in Ref. [93] can be replaced 
by the quantities Z^ff and the form factor F(— m^), where 



Zlf;F{-ml)a^ 
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em 



'Z 4:71 

and V^P^ /Z ~ F^") /N. The relevant box dia gram form factors are 
3 
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+ \Vtd\^ [Fxho^xt, Xi) - Fxbox(0, Xi) - Fxbox(0, Xt) + -Fxbox(0, 0)]} 



r>lJ.eqq 
^RR 
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{S ^V*; x^^ y,; xt^ Vt) 



where Xt = ml/my^^ and yt = ml/my^^. 

Finally we note that the presence of non-unitary mixing in the light neutrino sector (due 
to the matrix S ~ M^M^^) also affects the standard muon decay width, (and thus the 



-'^^ Although the process can also be mediated at tree-level by neutral Higgs bosons, these particles have to be 
very heavy due to constraints from K^-T^ mixing. 
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determination of G_f), as well as the capture rate for muons on the nucleus, F^apt- Exphcitly, 
one has 

r^^rW(l-[55t]ee-[55t]^^) and T^apt =^ rfj^, (l - [55^]^^) , (A-25) 

where T^/^^ and T^^Jpt are the SM values and we have omitted terms of order 5"'. These ex- 
pressions occur in the denominators of the branching ratio formulae in Eq. (58), and since the 
numerators are in general proportional to 0{S^) the effect will be negligible; in our analysis 
we use the standard value = G|,m^/(1927r^). 



A. 3. Loop functions 

The relevant loop functions are 
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